Two-Beam Spin Noise Spectroscopy 
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We propose a method of two-beam spin noise spectroscopy to test the spin transport at equilibrium via analysis 
of correlations between time-shifted spin fluctuations at different space locations. This method allows one to 
determine the strength of spin-orbit interaction and spin relaxation time and separate spin noise of conducting 
electrons from the background noise of localized electrons. We formulate a theory of two-beam spin noise 
spectroscopy in semiconductor wires with Bychkov-Rashba spin-orbit interaction taking into account several 
possible spin relaxation channels and finite size of laser beams. Our theory predicts a peak shift with respect 
to the Larmor frequency to higher or lower frequencies depending on the strength of spin orbit interaction and 
distance between the beams. The two-beam spin noise spectroscopy could find applications in experimental 
studies of semiconductors, emergent materials and many other systems. 
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The ability to understand and control the electron spin 
polarization in semiconductors will have profound tech- 
nological impacts 1,2 . For this purpose, it is important to 
obtain information about interactions that spins of con- 
ducting electrons experience in semiconductors as these 
interactions can play a significant role in the processes 
of spin relaxation and control. Experimentally, the elec- 
tron spin polarization and its fluctuations (spin noise) 
in semiconductors are often detected using the optical 
Faraday or Kerr rotation technique (see, e.g., Refs. 3-5). 
The Faraday rotation technique relies on the fact that 
a linearly polarized light propagating through a material 
with a net magnetic moment is subject to a rotation of its 
plane of polarization. The amount of rotation induced is 
proportional to the internal local magnetization of mate- 
rial and thus can provide information about electron spin 
polarization and its fluctuations. 

The major advantage of spin noise spectroscopy 4-8 to 
all other probes of semiconductor spin dynamics lies (e.g., 
optical pump-probe measurements 9 ) in the fact that in 
principle no energy has to be dissipated in the sample, 
i.e., it exclusively yields the intrinsic, undisturbed spin 
dynamics 5 . In addition, the spin noise spectroscopy al- 
lows accumulation of a large statistics that smoothes out 
the statistical noise in the data. At the present time, 
experimentalists utilize a single beam setup measuring a 
noise correlation function at the same space location 4-7 . 

We suggest a two-beam spin noise spectroscopy measur- 
ing both time and space correlations of spin fluctuations. 
Fig. 1 presents a possible experimental setup involving 
two laser beams (that may or may not overlap) separated 
by a distance d. This setup captures the propagation of 
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spin fluctuations between different space locations pro- 
viding a unique information regarding spin transport at 
equilibrium. To the best of our knowledge, there are no 
other available experimental approaches to obtain such 
information. 

In our theory, we assume that at the initial moment of 
time t — (that can be arbitrarily selected) the vector of 
spin polarization density is given by a vector of continu- 
ous random variables S(r, 0) = £(r) such that (£;(r)) = 
and (£i( r )£j( r ')) = X6(r — r')5y, where (..) denotes av- 
eraging over different realizations, i.j = x,y,z, and A 
is a parameter describing the strength of spin fluctua- 
tions. Using statistical considerations 10 , one can find 
that A = n/4, where n is 2D electron density. For a 
given realization of initial spin polarization density, the 
spin polarization at t > can be written using a Green's 
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FIG. 1. Two-beam spin noise spectroscopy. Correlations be- 
tween time-shifted intrinsic spin fluctuations at different space 
locations provide information about spin transport at equilib- 
rium. Here, d is the displacement of Beam 2 with respect to 
Beam 1. 
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function Gij of spin diffusion equation 



Si(r,t)= G ij (T,t;r , ,0)S j (r',Q)dv' + Xi {r,t), (1) 



where A is the sample area, Xj(r, f) is a stochastic func- 
tion describing the spin fluctuations created in the system 
at t > so that Xi( r ,0) = 0, and dr' = dx'dy' (in 2D 
case). Next, we note that each beam in the optical setup 
shown in Fig. 1 averages the space distribution of the 
Faraday rotation angle 6(r,t) — K,S z (r,t) according to 



,(*) 



f 



I m (v)6(r,t)dx 



K 



I m (r)S z (r,t)dr, 



m 

k is a constant that couples z-component of spin polariza- 
tion density with a local value of Faraday rotation angle, 
Po is the integrated laser beam intensity (power), and 
7 m (r) is the space distribution of m-th beam intensity, 
where m = 1,2. 

In the case of two-beam spin noise spectroscopy (we 
assume that two beams with the same incident polariza- 
tion and intensity are shifted by d as shown in Fig. 1), 
the experimentally determined correlation function is 



R(t) = (^(0)^(4)). 



(3) 



Using Eqs. (l)-(3), independence of x(r,i) from the ini- 
tial spin density distribution £(r), and the expression 
(6( r Kj( r ')) = ^( r — r ')<% (introduced above Eq. (1)), 
we find the general equation determining the second or- 
der spin noise correlation function: 



R(t) 
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^ J J J(r)/(r' -d)G^(r ) t;r',0)drdr / . 
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The Fourier transform of R(t) with respect to t is the 
noise power spectrum 



S(u) 



R(t) cos(art)d*. 



(5) 



We emphasize that although Eq. (4) contains only zz 
component of the Green function, the latter (as a solu- 
tion of a system of spin diffusion equations) incorporates 
both transverse and longitudinal dynamics of spin polar- 
ization. 

Eqs. (4) and (5) are quite general and can be ap- 
plied to many different experimental situations. Here, 
we consider a specific experimental system - a wire with 
the Bychkov-Rashba 11 spin-orbit interaction (SOI) and 
further develop the theory of two-beam spin noise spec- 
troscopy for this particular case. We assume that the 
wire is placed in an in-plane magnetic field perpendicu- 
lar to the wire axis (to shift a peak in the noise power 
spectrum from zero frequency) and, in addition to the 
D'yakonov-Perel' (DP) spin relaxation channel, there are 
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FIG. 2. Spin noise power spectrum in one-dimensional wires 
for different values of the beams' separation d. The direction 
of the peak's displacement is determined by the direction of 
the displacement of the second beam with respect to the first 
one along the wire. These curves were obtained using the 
parameter values Rorj = 1.57T, 7 = 0.2/tdp, and ujl = 5/tdp- 



additional spin relaxation mechanisms such as the Elliot- 
Yaffet mechanism 12,13 and/or electron spin relaxation on 
nuclear spins 14 . Assuming a Gaussian distribution of 
the incident laser beam intensities along the ^-direction, 
namely, I(x) oc exp(— x 2 /(2Rq)), where Ro is the beam 
radius, and using the Green's function of one-dimensional 
spin diffusion equation, 



G zz (x 1 1; x , 0) — 



1 



(x-a:') 2 
4Dt 



7 * cos (rj(x — x 1 ) — uj L t) , 
(6) 



where D is the diffusion coefficient, r\ is the spin preces- 
sion angle per unit length (due to the Bychkov-Rashba 
SOI), 7 is the spin relaxation rate due to additional to 
the DP spin relaxation channels, and wl is the Larmor 
frequency, we find (with a help of Eq. (4)) the noise 
correlation function 
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Unfortunately, the Fourier transform of Eq. (7) can 
not be explicitly obtained. It is possible, however, to 
derive a closed analytical expression in a limiting case 
when the short times provide the main contribution to 
the Fourier transform of Eq. (7). Strictly speaking, this 
limiting case is realized if Ro\r]\ S> 1 and d/(R^\r]\) <C 1. 
Neglecting Dt compared to Rq in Eq. (7), we rewrite this 
equation as 



R(t) oc cos 



dr]D\ 
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(8) 



The Fourier transform of Eq. (8) gives the spin noise 
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FIG. 3. Position of the peak in the spin noise power spectrum 
(Fig. 2) as a function of beams' separation d. These curves 
were obtained using the parameter values 7 — 0.2 /tdp and 
ul = 5/r D p- 
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FIG. 4. Half width at half-maximum (HWHM) for the peak 
in the spin noise power spectrum (Fig. 2) as a function of 
beams' separation d. These curves were obtained using the 
parameter values 7 = 0.2/tdp and ujl = 5/tdp- 



power spectrum 



S(uj) oc 



1 



(9) 



where r^ 1 = 7 + T^p is the total spin relaxation rate, 
Tr>p = (rfD)^ 1 is the D'yakonov-Perel' spin relaxation 
time 15,16 , and 5 = drjD/(Rg) is the peak shift from its po- 
sition at d = 0. According to Eq. (9), the spin noise spec- 
trum shows a Lorentzian peak centered at w' = uj l + 6; 
the resonance width is determined by the combined (due 
to all possible spin relaxation processes) spin relaxation 
time t s . Importantly, the resonance is shifted from ujl by 
S = dr/D/Rg depending on the strength of spin-orbit in- 
teraction (through the parameter 77) and the beams' sep- 
aration d. We emphasize that this shift of the resonance 
is an entirely new result that so far was not anticipated in 
the literature. We understand the resonance shift as fol- 
lows. Depending on the direction of electron diffusion (in 
the positive or negative x direction), the electron spins 
observed by the first beam acquire a spin rotation in the 
clockwise or counterclockwise direction (due to the spin- 
orbit interaction) by the moment of time when some of 
them are tested by the second beam in its own region of 
space. This spin rotation angle is added or subtracted 
from the Larmor precession angle causing a positive or 
negative peak shift. 

Fig. 2 represents the spin noise power spectrum ob- 
tained as the Fourier transform of Eq. (7) for two val- 
ues of d. Clearly, the peaks are of Lorentzian shape and 
drj = 47r peak is shifted by <5 with respect to drj = peak. 
While the half width at half-maximum (HWHM) gives 
information about the possible spin relaxation channels, 
the value of the peak shift 6 allows finding 77, so that the 
strength of the DP spin relaxation and its role in overall 
spin relaxation can be determined. 

Interestingly, additional effects can be observed if we 
consider a non-equilibrium situation e.g. by applying a 



weak external electric field E along the wire. The dif- 
fusion part of the propagator in (6) is then modified to 
■s/A-nDt ex P (~( x ~ x ' ~ (J'Et) 2 / iDtj , where fx is the elec- 
tron mobility. By neglecting ~ E 2 terms in correspond- 
ing R(t) function, the integration over the beam intensi- 
ties gives us 
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(10) 

According to this result, the effect of the electric field 
is twofold. First, it renormalizes the effective correla- 
tion time t s because spins moving along the electric field 
have a better chance to travel the distance d between the 
two beams before they relax, than spins moving against 
the electric field. This effect should be observable as 
the change of the peak amplitude as a function of the 
electric field. Perhaps more importantly, the electric 
field renormalizes the position of the peak maximum: 
S — > S — r\[iE. This shift can be used as an indepen- 
dent probe of the spin orbit coupling 77. It emerges 
because the nonzero average drift velocity fiE, via the 
spin orbit coupling, creates on average an effective trans- 
verse Zeeman-like field acting on electron spins with the 
corresponding Larmor frequency ~ r/fiE. This effect is 
not specific for quasi- ID geometry. It should emerge in 
any material with a linear in momentum spin-orbit cou- 
pling, such as the Rashba and the Dresslhaus couplings 
in 2D electron systems and strained 3D semiconductors. 
For example, taking parameters from Rashba electron 
system 17 : /i = a xx /e, where a xx = 1.8 • 10 3 il _1 m _1 is 
the conductivity, and 77 = 1.5 • 10 -13 eV-m, we find that 
an electric field E ~ lOeV/cm shifts the spin noise power 
peak by ~ 100MHz, which should be clearly observable 
in GaAs. 

Finally, we would like to consider 5 and HWHM be- 
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yond the limiting case discussed above. For this purpose, 
we perform exact calculations of noise power spectrum 
(based on Eq. (7)) and extract S and HWHM from these 
calculations. Figs. 3, 4 present some of our results. We 
note that while the peak shift scales almost linearly in 
the range of selected parameter values, the peak width 
demonstrates a significant variation at the smallest values 
Roij. Consequently, in the latter case the spin relaxation 
time should be extracted using the Fourier transform of 
the exact spin noise correlation function. 

In conclusion, as the peak shift and width can be eas- 
ily determined from experimental measurements, we an- 
ticipate that the suggested approach of two-beam spin 
noise spectroscopy will be extremely useful in the area of 
semiconductors and emergent materials. Possible appli- 
cations of the two-beam spin noise spectroscopy include 
the possibility to separate contributions of conduction 
and localized electrons to the spin noise, to identify the 
role of different spin relaxation mechanisms contributing 
to the process of electron spin relaxation in a particular 
material, and to measure spin orbit coupling anisotropy. 
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